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We review the properties of the nonlinearly dispersive Navier-Stokes- 
alpha (NS— a) model of incompressible fluid turbulence - also called 
the viscous Camassa-Holm equations in the literature. We first re-derive 
the NS— a model by filtering the velocity of the fluid loop in Kelvin's 
circulation theorem for the Navier-Stokes equations. Then we show 
that this filtering causes the wavenumber spectrum of the translational 
kinetic energy for the NS— a model to roll off as A: -3 for ka > 1 in three 
dimensions, instead of continuing along the slower Kolmogorov scaling 
law, fc -5 / 3 , that it follows for ka < 1. This rolloff at higher wavenumbers 
shortens the inertial range for the NS— a model and thereby makes it 
more computable. We also explain how the NS— a model is related 
to large eddy simulation (LES) turbulence modeling and to the stress 
tensor for second-grade fluids. We close by surveying recent results in 
the literature for the NS— a model and its inviscid limit (the Euler— a 
model) . 
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1 Introduction 

The energy in a turbulent fluid flow cascades toward ever smaller scales until it 
reaches the dissipation scale, where it can be transformed into heat. This cas- 
cade - creating fluid motions at ever smaller scales - is a characteristic feature 
of turbulence. This feature is also the main difficulty in simulating turbu- 
lence numerically, because all numerical simulations will have finite resolution 
and will eventually be unable to keep up with the cascade all the way to the 
dissipation scale, especially for complex flows, e.g., near walls and interfaces. 

The effects of subgrid-scale fluid motions occurring below the available res- 
olution of numerical simulations must be modeled. One way of modeling these 
effects is simply to discard the energy that reaches such subgrid scales. This 
is clearly unacceptable, though, and many creative alternatives have been of- 
fered. A prominant example is the large eddy simulation (LES) approach, see, 
e.g., lU, @. The LES approach is based on applying a spatial filter to the 
Navier-Stokes equations. The reduction of flow complexity and information 
content achieved in the LES approach depends on the characteristics of the 
filter that one uses, its type and width. In particular, the LES approach intro- 
duces a length scale into the description of fluid dynamics, namely, the width 
of the filter used. Note that the LES approach is conceptually different from 
the Reynolds Averaged Navier-Stokes (or, RANS) approach, which is based 
on statistical arguments and exact ensemble averages, rather than spatial and 
temporal filtering. After filtering, however, just as in the RANS approach, 
one faces the classic turbulence closure problem: How to model the filtered- 
out subgrid scales in terms of the remaining resolved fields? In practice, this 
problem is compounded by the requirement that the solution be simulated 
numerically, thereby introducing further approximations. 

This paper begins by reviewing a modeling scheme - - called here the 
Navier-Stokes- alpha model, or NS— a model (also called the viscous Camassa- 
Holm equations in - || ) — that introduces an energy "penalty" inhibiting 
the creation of smaller and smaller excitations below a certain length scale 
(denoted alpha). This energy penalty results in a nonlinearly dispersive 
modification of the Navier-Stokes equations. The alpha-modification ap- 
pears in the nonlinearity, it depends on length scale and we emphasize that it 
is dispersive, not dissipative. We shall re-derive the modified equations in Sec- 
tion 2 from the viewpoint of Kelvin's circulation theorem. As we shall show in 
section 3, this modification causes the translational kinetic energy wavenum- 
ber spectrum of the NS— a model to roll off rapidly below the length scale 
alpha as k~ 3 in three dimensions, instead of continuing to follow the slower 
Kolmogorov scaling law, k~ 5 ^ 3 . This roll off shortens the inertial range of the 
NS— a model and thus makes it more computable. This is the main new result 
of the paper, so we shall comment now on its main implications. 
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Since the energy spectrum rolls off faster, the wavenumber k = n a at which 
viscous dissipation takes over in the NS— a model must be lower than for the 
original Navier-Stokes equations. Hence, for a given driving force and viscosity, 
the number of active degrees of freedom Ndof for the NS— a model must be 
smaller than for Navier-Stokes. In section 3 we give a heuristic estimate of the 
number Ndof and compare it with the rigorous estimate derived in |J for the 
fractal dimension Df rac of the global attractor for the NS— a model. Namely, 

D frac < (N do} f 2 and N dof = {Ln a f ~ ^Re 3/2 , (1.1) 

where L is the integral scale (or domain size), K a is the end of the NS— a 
inertial range and Re = L 4 ^ 3 el! 3 / v is the Reynolds number (with NS— a energy 
dissipation rate e a and viscosity v). The corresponding number of degrees of 
freedom for a Navier-Stokes flow with the same parameters is 

N™ ee (L/£ Ko ) 3 ~ Re 9 /' , (1.2) 

where £ko denotes the Kolmogorov dissipation length scale. The implication 
of these estimates of degrees of freedom for numerical simulations that access 
a significant number of them using the NS— a model would be an increase in 
computational speed relative to Navier-Stokes of 

Thus, if a tends to a constant value, say L/100, when the Reynolds number 
increases - as found in - |J by comparing steady NS— a solutions with 
experimental data for turbulent flows in pipes and channels - then one could 
expect to obtain a substantial increase in computability by using the NS— a 
model at high Reynolds numbers. An early indication of the reliability of 
using these estimates to gain a relative increase in computational speed in 
direct numerical simulations (DNS) of homogeneous turbulence in a periodic 
domain is given in f?J. There, a computational speed up was achieved by using 
the NS— a model in DNS of turbulence in a periodic domain by a factor about 
equal to the fluctuation Reynolds number (« 250). In the case considered in 
0, this factor happens to be about Re/ 100. 

In section 4, we discuss the relation of the NS— a model to similar equa- 
tions which are derived in a different physical context, namely in the context of 
non-Newtonian fluids. The main difference between the NS— a model and the 
non- Newtonian second-grade fluids which it resembles lies in the choice of dissi- 
pation. The NS— a model uses the standard Navier-Stokes viscosity, while the 
second-grade fluid uses a weaker form of dissipation - namely, wavenumber in- 
dependent damping of the fluid velocity. Discussions of the relative advantages 
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of the two forms of dissipation (e.g., in terms of boundary data requirements 
and well-posedness) are beyond the scope of this review. However, we remark 
that the steady solutions of the NS— a model with the standard Navier-Stokes 
viscosity were found in 0] - to agree with experimental mean velocity pro- 
file data for turbulent flows in pipes and channels. The corresponding steady 
solutions of the second-grade fluid with its weaker form of dissipation were 
not found to so agree with this experimental data. Finally, in section 5, we 
provide a brief guide to the recent literature for those who might be interested 
in the mathematical context in which the alpha model was originally derived, 
as well as in its potential applications in turbulence modeling. In regard to 
the latter, the NS— a model was recently shown to be transformable to a gen- 
eralized similarity model for large eddy simulation (LES) turbulence modeling 



2 Kelvin-filtered turbulence models 

Although it was first derived from energetic considerations using the Euler- 
Poincare variational framework in M, [TTJ, the Navier-Stokes-alpha model may 
be motivated by an equivalent argument based on Kelvin's circulation theorem. 
The original Navier-Stokes (NS) equations are 

<9v 

— +v Vv + Vj9 = z/Av + f , with V • v = , (2.1) 
at 

for a forcing f and constant kinematic viscosity v. These equations satisfy 
Kelvin's circulation theorem, 

I v-dx= <f (z/Av + f ) • dx , (2.2) 

for a fluid loop 7(v) that moves with velocity v(x, t), the Eulerian fluid veloc- 
ity. 



Kelvin-filtering the Navier-Stokes equations. The equations for the 
NS— a model emerge from a modification of the Kelvin circulation theorem 
( |2.2| ) to integrate around a loop 7(11) that moves with a spatially filtered 
Eulerian fluid velocity given by u = g*v, where * denotes the convolution, 

u = g * v = J g(x - y) v d 3 y . (2.3) 

The "inverse" is denoted 



v = On, 



(2.4) 
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thereby defining an operator O whose Green's function is the filter g and which 
we shall assume is positive, symmetric, isotropic and time-independent. Under 
these assumptions the quantity (kinetic energy) 

E = i J u • v d 3 x = i J v • # * v d 3 x = ~ / u • Ou rf 3 x , (2.5) 

defines a norm. 

We obtain a modification to the Navier-Stokes equations (|2.1j) by replacing 
in their Kelvin's circulation theorem ( |2.2| ) the loop 7(v) with another loop 
7(11) moving with the spatially filtered velocity, u. Then we have, 

— I v ■ rfx = I (i/Av + f ) ■ dx . (2.6) 

/ 7 (u) / 7 (u) 

After taking the time derivative inside the Kelvin loop integral moving with 
filtered velocity u and reconstructing the gradient of pressure, we find the 
Kelvin- filtered Navier-Stokes equation, 

— + u ■ Vv + Vu T • v + Vp = z/Av + f , (2.7) 

with 



V-u = 0, and v = Ou . (2.8) 

The velocity u(x, t) is the spatially filtered Eulerian fluid velocity in equation 
( ^73|) . Note that continuity equation is now imposed as V • u = V ■ (g * v) = 0. 
The energy balance relation derived from the Navier-Stokes-alpha equations 
(|Jis 

j t yiu-vci 3 x = y u-f d 3 x- J u\V0 1/2 u\ 2 d 3 x, (2.9) 

where for the moment we have dropped the boundary terms that appear upon 
integrating by parts. That is, for the moment we ignore the boundary ef- 
fects and consider either the case of the whole space with solutions vanishing 
sufficiently rapidly at infinity, or the case of periodic boundary conditions. 

Similarities with previous work. 

1. Except for the term (Vu) ■ v, the Kelvin-filtered Navier-Stokes equation 
Q2.7|) is otherwise quite similar to Leray's regularization of the Navier- 
Stokes equations proposed in 1934 [|TTJ. Extension of the Leray regu- 
larization to satisfy the Kelvin circulation theorem was cited as an out- 
standing problem in Gallavotti's review |T!| . Looking at the equation of 
motion for the vorticity q = V x v, reveals even more similarity with 
Leray's regularization of the Navier-Stokes equations. (See the section 
about the vorticity below.) 
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2. At first glance, the Kelvin-filtered equations ( |2.7| ) in the absence of dissi- 
pation and forcing may seem reminiscent of a form of the Euler equations 
that was discussed in Kuz'min [13| and Oseledets (see also Gama 
and Frisch |1|), namely, 



dt dxi 



li 
duj 
dx i 



-li 

Ui + 

0. 



dv? 

dx i 
d<p 

dx % 



(2.10) 



In light of the third equation in this set, the second one is essentially 
the Hodge decomposition of the vector 7 with Euclidean components 
Comparing these equations with Euler's equations for the incompressible 
motion of an ideal fluid, 



dui • dui 

— + u>— + V P 







duj 
dx { 







(2.11) 



gives a relation between the pressure gradient and the "gauge function" 
4> that is reminiscent of (but different from) Bernoulli's law, 



Vl]|u| 



dt 



3 

U — — : 

OX 3 J 



(2.12) 



This type of relationship also arises in Hamilton's principle when one 
uses Clebsch variables [|T^] , in which case is a Lagrange multiplier that 
enforces the continuity equation. In contrast, equations ( |2.7| ) in the ideal 
unforced case give the Kelvin- filtered Euler model, 



dt dxi 

Vi 

duj 
dx l 



dv? 
dx % 



dp 

dx % 



(2.13) 



0. 



The Kuz'min-Oseledets form of the Kelvin-filtered Euler model appears 
by replacing Ui — > Vi only in the second of the three equations in (|2.10| ). 



Conservation laws. References to the interesting geometrical properties 
of the Euler— a model equations - namely equations ( |2.13| ) when O is the 
Helmholtz operator - are cited in the last Section. At this point, we only 
comment that in the case of periodic boundary conditions, or the case of 
the whole space with solutions vanishing sufficiently rapidly at infinity, these 
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nondissipative equations preserve the following two quadratic invariants, the 
kinetic energy, 



E= - I u-vd 3 x, (2.1 I) 



A = / v • curl vd 3 x. (2.15) 



and the v— helicity, 



The kinetic energy, E a , and the v-helicity, A, defined above, are also pre- 
served in bounded domains provided appropriate boundary conditions are im- 
posed. Boundary conditions sufficient for energy conservation when O is the 
Helmholtz operator 1 — a 2 A are 

iix(n-(Vu + Vu r ))=0, (2.16) 

where the vector n is normal to the boundary and superscript (-) T denotes 
matrix transpose. Boundary conditions sufficient for v-helicity conservation 
in general are A nu = and ph- curl v = 0, as seen from the helicity equation, 

d\ 

— = — div (A u + pcurl v) , where A = v-curlv. (2-17) 
at 

This equation is obtained by using only the motion equation in fl2.13|) and its 
curl. Therefore, conservation of helicity holds with these boundary conditions 
simply because of the form of the Kelvin-filtered motion equation in ( 2.131) , 



independently of the relation between v and u, and regardless of whether u is 
incompressible. 



Vortex transport and stretching. Let q = V x v be the vorticity of 
the unfiltered Eulerian velocity. The curl of the Kelvin-filtered Navier-Stokes 
equation ( |2.7| ) gives the vortex transport and stretching equation, 

^ + u- Vq-q- Vu = z/Aq + V x f . (2.18) 

at 

We note that the coefficient Vu in the vortex stretching term is the gradient 
of the spatially filtered Eulerian velocity u. Thus, Kelvin-filtering tem- 
pers the vortex stretching in the modified Navier-Stokes equations ( |2.7| ), while 
preserving the original form of the vortex dynamics. This tempered vortic- 
ity stretching is also reminiscent of Leray's approach [jnj] for regularizing the 
Navier-Stokes equations. 

Indeed, when O is the Helmholtz operator we take full advantage of this 
regularized vorticity stretching effect to prove in || the global existence and 
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uniqueness of strong solutions for the three dimensional NS— a model. In fact, 
when the forcing is in a certain Gevrey class of regularity (real analytic), then 
the NS— a solutions are also Gevrey regular, so that the Fourier coefficients of 
the solution decay exponentially fast, with respect to the wavenumbers, at a 
rate that increases as a 2 increases. See |[L7|, [pl| for the corresponding analysis 



of the Gevrey properties of the Navier-Stokes equations. This result suggests 
that the filtering, or smoothing of the velocity u due to the presence of alpha 
in the NS— a model could enhance the decay of its wavenumber spectrum 
and produce an earlier observation of this exponentially falling tail when the 
forcing is sufficiently smooth, even when the viscosity is small. 

Specializing to the Navier-Stokes-alpha model. The special case of the 
NS— a model emerges from the Kelvin-filtered Navier-Stokes equation (|2.7|) 
when we choose the operator O to be the Helmholtz operator, thereby intro- 
ducing a constant a that has dimensions of length, 

= l-a 2 A, with a = const. (2.19) 



In this case, the filtered and unfiltered fluid velocities in equation ( |2.4j ) are 
related by 

v = (l-a 2 A)u. (2.20) 

The Navier-Stokes-alpha model is given by the Kelvin-filtered Navier- 
Stokes equation ( p.7|) with definition Q2.20| ). The original derivation of the ideal 
Euler- alpha model (the v = case of the NS— a model) obtained by using 



the Euler- Poincare approach is given in ||, flO|| . The physical interpretations 
of u and v as the Eulerian and Lagrangian mean velocities are given in W§ . 
The corresponding kinetic energy norm ( |2.5| ) for the NS— a model is given 

by 



E, 



a 



1, 19 « 2 ,^ ,9 

- u 2 + — Vu 2 

2 2 1 



d 3 x, (2.21) 



and we repeat that u is the spatially filtered Eulerian fluid velocity. This 
kinetic energy is the sum of a translational kinetic energy based on the filtered 
velocity u, and a gradient- velocity kinetic energy, multiplied by a 2 . Thus, 
by showing the global boundedness, in time, of the kinetic energy ( |2.21|) one 
concludes that the coefficient Vu in the vortex stretching relation ( |2.18p is 



bounded in L 2 for the NS— a model. The second term in the kinetic energy 
( ^.21j ) imposes an energy penalty for creating small scales. The spatial integral 
of | Vu| 2 in the second term has the same dimensions as filtered enstrophy (the 
spatial integral of |V x u| 2 , the squared filtered vorticity). For a domain with 
boundary the spatial integral of | Vu| 2 in the second term in the kinetic energy 
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norm ( 2.21 ) is replaced by the integral of trace (4D • D) , where D is the strain 
rate tensor, D = (l/2)(Vu + Vu T ) in Euclidean coordinates p0 |. For the case 
we treat here, in Euclidean coordinates and in the absence of boundaries, all 
of these norms are equivalent for an incompressible flow. 



3 Spectral scaling for the NS— a model 

A preliminary scaling argument. Scaling ideas originally due to Kraich- 
nan |21]] for the case of two-dimensional turbulence may be applied to estimate 
the effect of the second term in the NS— a model energy ([2.21|) on the energy 
spectrum in the present case. For sufficiently small wavenumbers (ka <C 1, 
but fcL > 1) the first term in the energy E a in ( [2.21| ) dominates and the 
second term may be neglected. In this wavenumber region, the standard Kol- 
mogorov scaling argument for turbulence gives a k~ 5 ^ 3 spectrum by the usual 
dimensional argument for the inertial range, 

k=\L}-\ e a =\L] 2 \T}-\ E a (k) = \Lf[T}- 2 = e a a k\i{ ka < 1 , (3.1) 

so that a = 2/3 and b = —5/3. Conversely, for sufficiently large wavenumbers 
(ka ^> 1) and the second term dominates in E a . A preliminary scaling ar- 
gument would then indicate a k~ 3 spectrum for E a in the large wavenumber 
region, according to 

k = [L]-\ Va = [T]- 3 , E a (k) = [L} 3 [T]- 2 = v a a k b , if fca»l, (3.2) 

so that a = 2/3 and b = —3 in this case. Here r] a is the rate of dissipation of the 
J |Vu| 2 part of the kinetic energy (with the same dimensions as enstrophy). 
Therefore, in the wavenumber region near ka ~ 1 the kinetic energy spectrum 
may be expected to have a break in slope and roll off from k~ 5 ^ 3 to k~ 3 scaling 



The expectation from this preliminary scaling argument seems to be con- 
firmed in direct numerical simulations of the NS— a model, as shown in Figure 
1 below taken from |7j . Figure 1 shows the energy spectra resulting from three 
DNS of the NS— a model with mesh sizes of 256 3 for three cases: with a/L = 
(the Navier-Stokes equations), 1/32 and 1/8 for the same viscosity v = 0.001. 
The corresponding Taylor microscale Reynolds numbers R\ are reported as 
147, 182 and 279, respectively. The higher a (higher Reynolds number) flows 
are found numerically to have more compact energy spectra. 

A more refined argument . A more refined argument for the spectral scal- 
ing of the NS— a model shall now be given that depends explicitly upon prop- 
erties of the nonlinearity. Following [ffiH , let u K denote the component of u 



formed by the Fourier modes of fluid velocity u with wavenumbers in [k, 2k) 
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Figure 1: The DNS energy spectrum, E(k) = E a {k), versus the wave number k for three 
cases with the same viscosities, same forcings and mesh sizes of 256 3 for a — (solid line), 
1/32 (dotted line) and 1/8 (dotted-dash line). In the inertial range (k < 20), a power 
spectrum with /c~ 5 / 3 can be identified. For finite a, this behavior is seen to roll off to a 
steeper spectrum for k > 1/a. 

and similarly let v K denote the corresponding component of v, where k is 
well beyond the active wavenumbers of the driving force. The NS— a energy 
balance for this component is then 

\~tA u ^ v k) + v(-Au K , v K ) — T K T 2k , (3.3) 
2 at 

where T K and T 2k denote, respectively, the energy transfer rate at k from the 
low wavenumber modes w< to the high wavenumber modes u K + u > , and at 2n 
from -u< + U2 K to u > , where m< and w> are defined as 

u < = ^^Uj , and u> = Uj . (3.4) 

j<K j>2K 

In particular, 



-(B(u<, u<), u K ) + (B(u K + u>, v K + v>), u<) , 



(3.5) 
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where the bilinear operator B(u,v) is given by 

B(u, v) = -P„ (u x ( V x vfj , (3.6) 

and P a is the L 2 orthogonal projection (Leray projection), see, e.g., || for 
more mathematical details. For equation (|3.3|) to hold, it is essential that 
B(u,v) has the property, 

(B(u,v),u) = 0. (3.7) 

Note: property (|3.7|) does not hold for the enstrophy transfer rate of the NS— a 
model, so the NS— a model does not possess a Kraichnan type inertial range 
due to an enstrophy cascade. Time-averaging equation ( |3.3|) gives 



u((-Au K ,v K )} = (T K ) - (T 2k ). (3.8) 



Consequently, upon introducing the energy spectrum E a {n) , this time-averaged 
energy transfer equation implies 

VK 3 E a (K)~v k 2 E q (k)cIk~ (T K )- (T 2K ). (3.9) 

J K 

Hence, as long as the energy dissipation rate is small compared to the energy 
transfer rate, i.e., provided 

VK Z E a {n) < (T K ), (3.10) 
we have the inertial range condition, 

(T K ) ~ (T 2K ). (3.11) 
The total energy dissipation rate is estimated as 



v 



u ■ (-Av) d 3 x ) . (3.12) 

[0,LP 



Kraichnan JHJ] posits the following mechanism for the turbulent cascade: In 
the inertial range the eddies u K transfer their energy to the eddies u 2k in the 
time r K it takes to travel their length ~ 1/k. Their average velocity being 

/If „ \ 1/2 / f 2K EMdK \ 1/2 ( kEJk) \ 1/2 , n 

the eddy energy exchange, or turnover time t k will be 

r K ~ 1/{kU k ) . (3.14) 
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Consequently, the total energy dissipation rate is related to the k spectral 



energy density by using (|3.13|) as 



e« ~ r« 



r-iK pzk 5/2 

r- 1 / E a (k)dk~ K U K E a (k)dk~ 2 2U/2 £ ft (^ /2 , (3-15) 

which yields the following spectral scaling law for the NS— a inertial range, 



£«(«) 



el /3 (l + a 2 K 2 f' 3 



K 



5/3 



(3.16) 



Thus, the total energy present in the NS— a inertial range is actually enhanced 
by the presence of alpha. Notice, however, that in terms of the filtered velocity 
u alone, the translational kinetic energy spectrum is given by 



EJk) 



2/3 



2/2, 

if KOt <ti 1 , 



K 



[1 + a 2 K 2 ) k 5 /3 (i + a 2 K 2 ) 2 / 3 



2/3 



(3.17) 



if na ^> 1 . 



Hence, the anticipated k~ 3 behavior appears in the spectrum for the a— filtered 
velocity as a consequence of equation p. 17] ) arising from Kraichnan's argument 
|2l| associating energy transfer rates and eddy turnover times for the filtered 
velocity. According to this argument, the presence of alpha reduces the energy 
associated with the higher wave numbers in the L 2 norm of the filtered veloc- 
ity u. This conclusion from (ft.lTp agrees with the trends shown in Figure 1 
obtained from high resolution DNS studies of the NS— a equation in [0. How- 
ever, these numerical studies do not have sufficient dynamic range to confirm 
this conclusion entirely. Thus, uncertainty remains in claiming a numerical 
confirmation because the quantity E a (n)/(1 + a 2 K 2 ) in equation (|3.17|) is un- 
affected by the alpha-modification for na 1 and may already be out of the 
inertial range and into the dissipation range for Ka ^> 1. Studies in progress 
using DNS of the limit a — >• oo NS— a equation hope to clarify this point 
Next we shall discuss the extent of the inertial range for the NS— a model. 
Equation (|3.17| ) holds only in the inertial range, that is, provided, cf. (|3.10| ), 



VK 3 E a (n) < e c 



(3.18) 



This inert ial-range inequality may be expressed equivalently as k C K a where 
n a (the end of the NS— a inertial range) is given in terms of the NS— a Kol- 
mogorov wavenumber n a ,Ko by using ( |3.15| ) and fl3.17| ) to find 



1 + a K a ) ~ -j 



h a,Ko 



(3.19) 
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Thus, as expected from the numerical simulations of |7J , the inertial range 
is shortened to k < n a for the NS— a model by its nonlinear dispersive filter- 
ing with lengthscale a. For sufficiently large NS— a Kolmogorov wavenumber 
Ka,Ko and with a fixed, the wavenumber n a at the end of the NS— a inertial 
range is determined from ( |3.19| ) to be 



1/3 



*> - {-) *Tko ■ (3-20) 

This is a relationship among the three progressively larger wavenumbers, 

1/a < n a < K a>Ko . 

Shortening the inertial range for the NS— a model to k < K a rather than 
k < K a< Ko implies fewer active degrees of freedom in the solution. 

Counting degrees of freedom. If one expects turbulence to be "extensive" 
in the thermodynamic sense, then one may expect that the number of "active 
degrees of freedom" N do j for alpha-model turbulence should scale as 

N dof ~ (Ln a ) 3 ~ (L/a)(LK a>Ko ) 2 ~ ^Re 3 ' 2 , (3.21) 

where L is the integral scale (or domain size), K a is the end of the NS— a 
inertial range and Re = L 4 ^ 3 eH 3 /u is the Reynolds number (with dissipation 
rate e a and viscosity v). The corresponding number of degrees of freedom for 
Navier-Stokes with the same parameters is 

N$ = (L K a>Ko f ~ Re 9 /\ (3.22) 

and one sees a possible trade-off in the relative Reynolds number scaling of the 
two models. Should these estimates of the number of degrees of freedom needed 
for numerical simulations using the NS— a model relative to Navier-Stokes not 
be overly optimistic, the implication would be one factor of (N^f j Ndof) 1 ^ 3 in 
relative increased computational speed gained by the NS— a model for each 
spatial dimension and yet another factor for the accompanying lessened CFL 
time step restriction. Altogether, this would be a gain in speed of 



4/3 / „, \ 4 / 3 

Re. (3.23) 



\N dof ) \L) 

Since a/L < 1 and Re 3> 1, the two factors in the last expression compete, 
but the Reynolds number should eventually win out, because Re can keep 
increasing while the number oc/L is expected to tend to a constant value, 
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say a/L = 1/100, at high (but experimentally attainable) Reynolds numbers, 
at least for simple flow geometries. Empirical indications for this tendency 
were found in - |J by comparing steady NS— a solutions with experimental 
mean- velocity-profile data for turbulent flows in pipes and channels. 

Thus, according to this scaling argument, a factor of 10 4 in increased com- 
putational speed for resolved scales greater than a could occur by using the 
NS— a model at the Reynolds number for which the ratio n a ,Ko/K a = 10. 
An early indication of the feasibility of obtaining such factors in increased 
computational speed was realized in the direct numerical simulations of homo- 
geneous turbulence reported in 0, in which n a ,Ko/ K a — 4 and the full factor 
of 4 4 = 256 in computational speed was obtained using spectral methods in a 
periodic domain at little or no cost of accuracy in the statistics of the resolved 
scales. 



Related mathematical results. The paper |J shows that strong solutions 
of the NS— a model exist globally, they are unique, and they lie on a global 
attractor whose fractal (Lyapunov) dimension is bounded above by 

V frac < NH 2 f , (3.24) 



with Ndof defined as in equation ( |3.21 ). This rigorous mathematical bound 



exceeds the expected value obtained above by heuristically counting degrees 
of freedom assuming "extensive" turbulence. Thus, it would imply a smaller 
increase in computational speed than that of Q3.23|) . However, this rigorous 
bound may have room for improvement. 



Oboukov cascade rate. We have established in equation ( 3.14| ) that the 



eddy turnover rate based on the translational velocity for the NS— a model is 
r K ~ kU k ~ \l ~ ( ^—^3 ) • ii-2->) 

Therefore, 



1 + a 2 K 2 \ 1 + a 2 K 2 



for low wavenumbers (not 1), information (or error) propagates be- 
tween scales K and 2k at a rate proportional to k 2//3 , in agreement with 



the classical accelerated cascade of Oboukov |25[ as cited, e.g., in |[26| ; 
while 

for high wavenumbers (kq ^> 1), information propagates at a constant 
rate, independently of wavenumber for the NS— a model (as occurs 
also in the case of 2D turbulence discussed by Leith and Kraichnan f27|). 
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Without an accelerated cascade at high wavenumber, the NS— a model should 
tend to be more predictable than the original Navier-Stokes model. Of 
course, it makes sense that an averaged, or filtered description of fluid flow 
would propagate high wavenumber information and errors at a slower rate 
than an instantaneous, unfiltered description does (e.g., think of climate vs 
weather). The n~ 3 spectrum of the translational kinetic energy in the alpha 
models for an ^> 1 is consistent with this interpretation of reduced error 
propagation rate as measured in the L 2 norm of the filtered velocity. 



Rheology of NS— a turbulence: 2nd-grade 
fluids 



We rewrite the NS-alpha equations (|2.7|) with v 
in their equivalent constitutive form, 



(1 — a 2 A)u and a constant 



du 

~dt 



V ■ T , where T = —pi + 2i/(l - a 2 A)D + 2a 2 D , 



(4.1) 



with V ■ u = 0, strain rate D = (1/2) (Vu + Vu T ), vorticity tensor Q = 

o 

(l/2)(Vu — Vu ), and co- rotational (Jaumann) derivative given by D = 
dD/dt + Dfi - VLD. 

In equation ( |4.1|) , one recognizes the constitutive relation for the NS— a 
model as a variant of the rate-dependent incompressible homogeneous fluid of 
second grade |2B |, p9fl , in which the dissipation, however, is modified by com- 
position with the Helmholtz operator (1 — a 2 A). Thus, the NS— a model has 
nearly the same stress tensor as the second-grade fluid. However, it is not quite 
the same. The stress tensor for the NS— a model uses Navier-Stokes viscous 
dissipation, instead of the weaker form of dissipation used for second-grade flu- 
ids that is independent of wavenumber. Note: despite first appearances, there 
is no hyperviscosity in the NS— a model, only the standard Navier-Stokes vis- 
cosity. 

Equations for the second grade fluid were treated recently in the mathemat- 
ical literature |3l| , pT 



32 



Also, in |M| local well-posedness results for some 
initial value problems for differential type fluids arrising in non-Newtonian 
fluids (second-grade and third-grade fluid) were obtained by transferring the 
problem from the Eulerian to the Lagrangian setting, thereby extending the 
Arnold [33] and Ebin-Marsden program [M\ to the case of second-grade flu- 
ids. For the case of second-grade fluids, for ka > 1 the wavenumber spectrum 
should still roll over to k~ 3 , provided the Kolmogorov/Oboukov argument still 
holds for the weaker dissipation. 

The association of turbulence closure models with non-Newtonian fluids 



is natural. There is a tradition at least since Rivlin |55] of modeling turbu 
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lence by using continuum mechanics principles such as objectivity and material 
frame indifference (see also f36|). For example, this sort of approach is taken 
in deriving Reynolds stress algebraic equation models [PTfl . Rate-dependent 
closure models of mean turbulence have also been obtained by the two-scale 
DIA approach |38[] and by the renormalization group method |39[ . 

Despite its similarity to the rheology of second-grade fluids, the alpha pa- 
rameter in the NS— a model is actually not a material parameter. Rather, it is 
a flow regime parameter. Since the NS— a model describes mean quanti- 
ties, it was proposed as a turbulence closure model and this ansatz was tested 
by comparing its steady solutions (using the standard Navier-Stokes viscosity) 
to mean velocity measurements in turbulent channel and pipe flows in ||-||. 
These experimental tests show that alpha depends slightly on Reynolds num- 
ber and varies slightly with distance from the wall, for the low to moderate 
Reynolds numbers available in channel flow. However, at the high to very 
high Reynolds numbers available in pipe flows, alpha becomes independent 
of Reynolds number and takes a small value - about one percent of the pipe 
diameter. 



5 Guide to recent NS— a model literature 

In f3j-p the authors introduce random fluctuations into the description of the 
fluid parcel trajectories in the Lagrangian in Hamilton's principle for ideal in- 
compressible fluid dynamics. They then take its statistical average and use the 
Euler-Poincare theory of ||, [|T(J], ]|0| to derive Eulerian closure equations for 
the corresponding averaged ideal fluid motions. The Euler-Poincare equation 
that is used in deriving the ideal dynamics of the NS— a model is equivalent 
in the Eulerian picture to the corresponding Euler-Lagrange equation for fluid 
parcel trajectories for Lagrangians that are invariant under the right-action of 



the diffeomorphism group. See |J, JjlJ], [2(J, and references therein for more 
discussions of Euler-Poincare equations. The Euler-Poincare theory is applied 
for modeling fluctuation effects on 3D Lagrangian mean and Eulerian mean 



fluid motion in 19 



The NS— a model equations (also sometimes called in the literature the 
viscous Camassa-Holm equations) are proposed in as a turbulence clo- 

sure approximation for the Navier-Stokes equations. The analytic form of the 
velocity profiles based on the steady NS— a equations away from the viscous 
sublayer, but covering at least 95% of the channel, depends on two free param- 
eters: the flux Reynolds number R, and the wall-stress Reynolds number R Q . 
(Due to measurement limitations, most experimental data are contained in 
this region.) The authors further reduce the parameter dependence to a single 
free parameter by assuming a certain drag law for the wall friction D ~ R^fR 2 . 
For most of the channel the steady NS— a solution is shown to be compatible 
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with empirical and numerical velocity profiles in this subregion. The NS— a 
steady velocity profiles agree well outside the viscous sublayer with data ob- 
tained from mean velocity measurements and simulations of turbulent channel 
and pipe flow over a wide range of Reynolds numbers, H-||. 



DNS and comparisons with LES. In |7j direct numerical simulations 
(DNS) of the NS— a equations are compared with Navier-Stokes DNS and 
interpreted as behaving like an LES (Large Eddy Simulation) model. That is, 
the NS— a model is shown to produce an accurate dynamical description of the 
large scale features of turbulence driven at large scales, even at resolutions for 
which the fine scales are not resolved. The NS— a-model may at first appear 
not to be an LES model, since it has rate-dependence that is not admitted 
by LES models. However, the NS— a model was recently shown to transform 
to a generalized LES similarity model in ||. This is a promising result for 
LES modeling, since the mathematical theorems available in || for existence, 
uniqueness and finite dimensional global attractor for the NS— a model will 
now be transferable to the continuous formulations of this class of generalized 
LES similarity models. As a basis for numerical schemes, the NS— a model 



also resembles vorticity methods, as observed in [PD|1 , pT| , 



Geodesic motion. The completely integrable one-dimensional Camassa- 



Holm equation [43] is expressed on the real line as 



oc 



^ + u^ + 2v^ = 0, «(*,*) = !/ e^v(y,t)dy. (5.1) 



Thus, we have v = u — d 2 u/dx 2 , cf. equations ( |2.7| ) with definition ( |2.20| ). This 
one-dimensional equation is formally the Euler-Poincare equation for geodesic 
motion on the diffeomorphism group with respect to the metric given by the 
mean kinetic energy Lagrangian, which is right-invariant under the action of 
the diffeomorphism group. See ||, [TIJ for detailed discussions, applications 



and references to Euler-Poincare equations of this type for ideal fluids and 
plasmas. See ||, [M] for the original derivation of the n-dimensional Camassa- 



Holm, or Euler— a equation in Euclidean space. See jPJ, [[45| for discussions of 
its generalization to Riemannian manifolds, its existence and uniqueness on a 
finite time interval, and more about its relation to the theory of second grade 
fluids. Additional properties of the Euler— a equations, such as smoothness of 
the geodesic spray (the Ebin-Marsden theorem) are also shown to hold in ME\ 



and the limit of zero viscosity for the corresponding viscous equations is shown 
to be a regular limit, even in the presence of boundaries for either homoge- 
neous (Dirichlet) boundary conditions, or for boundary conditions involving 
the second fundamental form of the boundary. Functional-analytic studies of 
the ideal Euler— a model are made in PO], ||46|| . The methods introduced and 



NS-alpha turbulence 



Foias, Holm & Titi 



18 



applied in j2(J, j46|, while geometrical in nature, also address analytical issues 
and obtain results that were not previously available by traditional techniques 
for partial differential equations. For example, these methods produce local in 
time existence of C°° viscosity independent solutions for the Euler— a equa- 
tions in n dimensions (and in particular 3D) for a fluid container that can be 



an arbitrary Riemannian manifold with boundary fi? 



Mathematical estimates for strong solutions with dissipation. Pa- 
per || provides the mathematical estimates that are needed to show that the 
solutions of the NS-ct model exist globally, are unique and possess a global at- 
tractor with finite fractal dimension satisfying equation Q3.24 ) . The estimates 
in do depend on the viscosity remaining positive. 
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